Ferroelectricity in the Dion-Jacobson CsBiNb207 from first principles 
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We have studied ferroelectricity in Dion-Jacobson CsBiNb207 from first principles. Using group- 
theoretical analysis and first-principles density functional calculations of the total energy and 
phonons, we perform a systematic study of the energy surface around a paraelectric prototypic 
phase. Our results suggest that CsBiNb207 is a ferroelectric with a polarization of Pa=40/iC cm"'^. 
We propose further experiments to clarify this point. 

PACS numbers: 77.80.Bh, 61.50.Ks, 63.20.Dj 
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The ABO3 perovskites with their corner-shared oc- 
tahedra, display an amazing variety of physical prop- 
erties by simply changing the A and B cations. In 
the rational design of new materials much effort has 
been directed towards manipulating the connectivity 
of perovskites by forming structures consisting of 2- 
d layers of perovskites separated by cation network. 
Homologous series of these layered perovskites such 
as the Dion-Jacobson A[A'„_iB„03„+i], Ruddlesden- 
Popper A2[A'„_iB„03„+i], and Aurivillius [Bi202] [ 
A'„_iB„03„+i] (where A=alkali metal, A' alkaline 
earth, and B transition metal) have been formed by a 
number of novel growth techniques including molecular 
beam epitaxy-'^ and soft chemistry techniques'^ such as 
ion-exchange and intercalation reactions. 

While ferroelectricity is ubiquitous in the ABO3 per- 
ovskite titanates, niobates, and tantalates, the only lay- 
ered perovskites to display a ferroelectric (FE) transi- 
tion has been the Aurivillius compounds, while a pre- 
dicted Ruddlesden-Popper FEJ^ has yet to be confirmed 
experimentally. Recently, Lightfoot and coworkers per- 
formed neutron diffraction on ceramic samples of the 
n—2 Dion-Jacobson compound CsBiNb2 07 and found 
that the room temperature structure crystallizes in the 
polar space group P2iami4 Nevertheless, from electrical 
measurements of the static dielectric constant, they con- 
cluded that CsBiNb207 does not display ferroelectricity. 

In this letter we perform group theoretical analysis and 
first-principles calculations of the structural energetics 
and phonons of CsBiNb2 07. We show that the polar 
structure observed in Ref. U is a stable low-symmetry 
phase of CsBiNb207 for which we calculate a substantial 
polarization, Pg^iOfiC cm~^. Further we point out that 
the observed polar distortions from a non-polar reference 
structure are consistent with the criteria established by 
Abrahams^ for systems having a high probability of dis- 
playing a FE transition. All this strongly suggests that 
CsBiNb207 is a FE. 

First-principles density-functional calculations using 
projector augmented-wave potentials were performed 
within the local density approximation (LDA) as im- 
plemented in the Vienna ah initio Simulation Package 
(VASP)i^ The wavefunctions were expanded in plane 
waves up to a kinetic energy cutoff of 500 eV. Integrals 
over the Brillouin zone were approximated by sums on 



a 6 X 6 X 2 F-centered /c-point mesh. Polarization was 
calculated using the modern theory of polarization^ as 
implemented in VASP. 

We approach the problem of identifying stable, low- 
symmetry phases of CsBiNb207 and possible ferroelec- 
tricity by investigating the lattice instabilities of a hy- 
pothetical reference structure. Reference structures of 
this type can be chosen by detecting a pseudosymmetrjiii 
of the observed polar phase, that is by finding structures 
in a group-subgroup relation where the atomic distortions 
relating the two phases are relatively small. The simple 
tetragonal (paraelectric) space group P4/mmm, which 
the related compound, CsLaNb207 has been shown to 
form at room temperaturai? satisfies this requirement. 
While it is reasonable to assume that at a high temper- 
ature CsBiNb207 transforms to P4/mmm, we must em- 
phasize that regardless of whether this transition occurs, 
P4/mmm is a good reference structure. 

From x-ray and neutron powder diffraction, Ref. 
showed that the unit cell of CsBiNb207 was orthorhom- 
bic at room temperature and in a \/2at x \/2at x ct rela- 
tion to the simple tetragonal reference structure, system- 
atic reflections being consistent with the non-polar space 
group Pmam and with the polar space group P2iam. 
Subsequent Rietveld refinements gave a better fit to the 
P2iam space group. 

Both Pmam and P2iam are subgroups of P4/mmm. 




FIG. 1: Structure of ferroelectric CsBiNb207 
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FIG. 2: Unstable phonons of the prototypic phase, P4/mmm, 
connecting possible low-temperature paraelectric (PE) and 
ferroelectric (FE) phases. 



TABLE I: Crystal structure of ferroelectric CsBiNb207, 
Space Group: P2iam, a= 5.412 A, 6=5.343 A, c= 11.22 A 
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Group-theoretic analysis shows that the zone-center F- 
phonons and the zone-boundary M-point phonons drive 
aU relevant transitions. Performing phonon calcu- 
lations at these selected points in the Brillouin zone 
and subsequently investigating the structural energet- 
ics around the reference structure facilitates a system- 
atic method of determining possible low-symmetry struc- 
tures. Specifically, as seen from Fig. |21 Pmam is re- 
lated to P4/mmm by the freezing-in of a single zone- 
boundary phonon. On the other hand, a F^ zone-center 
ferroelectric mode and a M-point zone-boundary mode 
are required to account for the symmetry change associ- 
ated with P2iam. 

We performed full optimization of the lattice parame- 
ters and internal coordinates of CsBiNb207 in the tetrag- 
onal reference structure space group P4/mmm. Our cal- 
culated structural parameters are displayed in Table ^ 
where they are written in terms of the V2at x V2at x c* 
lattice vectors so as to readily compare with the P2iam 
structure. Next we calculated the F^ phonons by com- 
puting the relevant block of the dynamical matrix from 
finite differences of the Hellmann-Feynman forces. We 
found one unstable Fg" phonon with frequency lo— i 185 
cm~^. The character of the real-space eigendisplacement 
of this unstable mode consists mostly of Bi-ions moving 
against oxygen in Bi-0 planes as would be expected due 
to the Bi 6s lone pairJ^ but also consists of a smaller 
amount of Nb-0 displacements. The freezing-in of this 
ferroelectric mode would lead directly to the C2mm space 
group without a doubling of the unit cell. Next we com- 
puted the M^, MJ, and phonons at the q=(l/2, 
1/2, 0) point of the Brillouin zone. has two unstable 
modes with frequencies ui— i 210 cm~^ and 102— i 70 
cm^^. The Mg" modes are antiferrodistortive, the char- 
acter of which are Bi-0 and Nb-0 displacements for uii 
and UJ2 respectively. The freezing-in of these Mg' modes 
leads to the paraelectric, doubled unit cell, Pmam space 
group. Finally, the and the phonons both con- 
sist only of oxygen motion in the Nb-0 planes. The 
phonon is unstable (w= i 190 cm^^) leading to space 
group P4/mbm while the phonon is stable. 

Based on the results of the phonon calculations, we see 
that transitions out of the P4/mmm reference structure 
to the space groups C2mm, Pmam, P4/mbm, and P2iam 



are all possible due to the freezing-in of one or more un- 
stable modes, as summarized in Fig.|5| To identify which 
phase is most energetically favorable and therefore most 
likely to be observed at low temperatures, wc performed 
a series of structural relaxations within each of the above 
mentioned space groups. We relaxed all lattice parame- 
ters and internal coordinates, the results are as follows: 



Space Group 
P4/mbm 
C2mm 
Pmam 
P2iam 



Irrep A E/ f.u. 
M+ -400 meV 
F.^ -550 meV 
-350 meV 



M+ 
and 



or T7 and 



-800 meV 



where it is clear that P2iam, being 800 meV lower than 
the tetragonal reference structure and 450 meV lower 
than Pmam, is the computed ground state, in full agree- 
ment with the Rietveld refinements of Ref. 0[ The re- 
sults of the structural relaxations within the P2iam space 
group are shown in Table |l| We find lattice constants 
«1% smaller than those determined experimentally at 
300 K, typical of LDA calculations. In Table ITU wc com- 
pare with the experimentally determined bond lengths. 
We find excellent agreement, within 1-2% except for the 
Bi-Oa bonds which differ by 3-5%. 

From phonon calculations, we have shown that lattice 
instabilities of a high symmetry reference structure of 
CsBiNb207 lead naturally to a P2iam ground state. As 
we have stated the symmetry of this space group allows a 
ferroelectric polarization to develop along the a-axis. In 
fact, from Berry-phase calculations of the polarization we 
find that CsBiNb207 should have a polarization P^w 40 
/iC cm~^, which is about three times larger than the po- 
larization of the n=2 Aurivillius compound SrBi2Ta209 
and comparable to that of the n—3 Bi4Ti30i2. 

In principle, a material is considered a ferroelectric if 
both 1) it displays a spontaneous polarization (P^) and 
2) if Ps can be flipped in an electric field. This defi- 
nition is reflected in a experimentally derived structural 
criteria for identifying possible ferroelectrics as 1) the 
largest Ax > 0.1 A and 2) all Ax < 1.0 A, where Ax 
can be though of as an atomic distortion along the polar 
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TABLE II: Bond lengths (A) of CsBiNbaOT, Space Group; 
P2iam, Exp.: a= 5.495A, 6=5.423A, c= 11.38A. LDA: a= 
5.412A, &=5.343A, c= II.22A 



Atoms Exp. Ref. 4° LDA 

Cs - 0(2) 3.038 2.968 

3.144 3.096 

3.202 3.133 

3.232 3.184 

Nb - 0(1) 2.403 2.375 

Nb- 0(2) 1.755 1.764 

Nb- 0(3) 1.941 1.957 

2.024 2.037 

Nb- 0(4) 1.973 1.948 

2.082 2.031 

Bi - 0(i) 2.257 2.224 

2.269 2.243 

3.243 3.199 

3.256 3.228 

Bi - 0(3) 2.707 2.614 

2.794 2.654 

Bi - 0(4) 2.409 2.398 

3.311 3.354 



"Powder neutron diffraction, room temp. 

axis from the centrosymmetric reference structured In 
CsBiNb207, the largest atomic distortion relative to the 
oxygen framework is Ax(Bi)R:i0.5A which clearly meets 
Abrahams' structural criteria and suggests, independent 
of the first-principles results, that CsBiNb207 is a FE. 

Experimentally, an indication of a ferroelectric transi- 
tion is usually shown by the strong temperature depen- 
dence of the static dielectric constant around Tc, nor- 



mally a Curie-Weiss behavior. At temperatures much 
lower than Tc, a strong temperature dependence of the 
dielectric constant is not expected, for example, the static 
dielectric constant of LiNbOa {Tc= 1450 K) over the tem- 
perature range 0K-300K. The fact that the static dielec- 
tric constant of CsBiNb2 07 does not display a large tem- 
perature variation around room temperature does not 
preclude a ferroelectric transition at a higher tempera- 
ture. Considering the PE-to-FE energy scale that we 
calculated, AE«-400meV/f.u., Tc is expected to be much 
higher than room temperature. This is consistent with 
other known ferroelectrics where the stereochemical ac- 
tivity of the Bi lone pair is the driving force towards fer- 
roelectricity such as in BiFeO^i^ and in SrBi2Ta209iii In 
fact, from an empirical expression that has proven quite 
successful in predicting Tc of several new ferroelectrics, 
we find Tc = 2 x W^iAx)^ K ^ 810 K where Ax is the 
largest metal atom (in this case Bi) displacement along 
the ferroelectric direction after a shift of origin so that 
Ax(Cs)-hAx(Bi)-hAx(Nb)=Oia We suggest that direct 
measurements of the polarization in CsBiNb207 be per- 
formed, such as E-P hysteresis loops, to reveal it's ferro- 
electric nature since the actual Tc may in fact be higher 
than the decomposition temperature of the material. 

We have confirmed that the ground state of 
CsBiNb207 forms in the polar space group P2iam and 
expect a substantial polarization P^ « 40 /xC cm~^. The 
similarity of the computed lattice instabilities to the n—2 
Aurivillius compoundsii is intriguing, the implications of 
which should be pursed further given the fundamental 
and technological interest in such materials. The Dion- 
Jacobson family of compounds may provide another av- 
enue towards the design of new ferroelectrics and, given 
the propensity towards ion-exchange, multiferroics^ 
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